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Abstract 

We consider a model for a flat, disk- like galaxy surrounded by a 
halo of dark matter, namely a Vlasov-Poisson type system with two 
particle species, the stars which are restricted to the galactic plane and 
the dark matter particles. These constituents interact only through 
the gravitational potential which stars and dark matter create collec- 
tively Using a variational approach we prove the existence of steady 
state solutions and their nonlinear stability under suitably restricted 
perturbations. 

1 Introduction 

Around 1970 astrophysicists noticed that in typical spiral galaxies the rota- 
tion velocities of the stars, when computed in the gravitational potential of 
the visible matter, do not fit with their observed rotation velocities. It was 
then conjectured that such galaxies are surrounded by a halo of so far not 
directly observable dark matter in such a way that the rotation velocities of 
the stars are consistent with the resulting gravitational potential [8]. For an 
introduction to dark matter we refer to [31 Chapter 10] and the references 
there. 

The distribution of the stars in a galaxy is usually modeled by a density 
function on phase space, and it is assumed that collisions are sufficiently rare 
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to be neglected and that the stars interact only by the gravitational potential 
which they create collectively. In a non-relativistic setting this results in a 
system of partial differential equations which in the mathematics literature is 
known as the Vlasov-Poisson system, cf . [22] . While the true physical nature 
(and existence) of dark matter are still conjectural, we are aware of at least 
one astrophysics investigation where it is also modeled as Vlasov-type matter, 
cf. [23]. Given the fact that the only role which galactic dark matter has to 
play is to provide the mass and hence the gravitational potential needed to 
resolve the discrepancy concerning the rotation velocities of the stars, such 
a description of dark matter seems natural. 

In the present paper we investigate a model for a flat, disk-like galaxy 
with a halo of dark matter where both the distribution of the stars in the 
galactic plane and the distribution of the dark matter particles in the halo 
obey a Vlasov equation, and the interaction among stars, dark matter, and 
between these two constituents is through the gravitational potential which 
all the particles (stars and dark matter) create collectively. 

Following the practice in astrophysics we assume that the stars are re- 
stricted to a plane which we take to be the Xi,x 2 plane. Their distribution on 
phase space is given by / = f(t,x,v) > where t > denotes time and x,v G 1R 2 
denote position and velocity in the galactic plane. The distribution of the 
dark matter particles is given by f = f(t,x,v) >0 where x,v GlR 3 denote po- 
sition and velocity in three dimensional space. The evolution of the galaxy 
and its halo is then governed by the following Vlasov-Poisson type system of 
equations: 

d t f + v-V x f-V x U e -V v f = 0, (1.1) 
dJ+v-Vj-V x U e (;0)-Vj = 1 (1.2) 

U e (t,x) = U(t,x) + U(t,x) = - [ p^rdy- [ , ^^ dy, (1.3) 

p(t,x) = / f(t,x,v)dv, p(t,x)= / f(t,x,v)dv. (1.4) 
Jm 3 Jr 2 

Here p and p are the spatial mass densities of dark matter respectively stars, 
U and U are the induced Newtonian potentials, and U e denotes the poten- 
tial of the system as a whole, i.e., the effective potential which determines 
the particle orbits. In order that the stars remain in their plane it is suffi- 
cient to require that f(t,x,X3,v,vs) —f(t,x, — xs,v,— V3), a condition which at 
least formally is preserved by solutions of the system and which implies that 
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VU(t,x,0) is parallel to the plane; for VU(t,x) this is true automatically. 
Throughout this paper we use the convention that variables with (without) 
tilde denote flat (non-flat) quantities. 

To our knowledge a fully non-linear model where the gravitational inter- 
action within both types of matter and between the two types is taken into 
account has so far not been investigated. Our aim is to prove the existence 
and non-linear stability of steady state configurations to this system. We 
obtain such stable steady states as minimizers of the total energy 

2 [I \v\ 2 fdxdv + - II \v\ 2 fdxdv 



If 1 f 

+- I U e (x)p(x)dx + — / U e (x,0)p(x)dx, 



satisfying suitable constraints. This so-called energy- Casimir approach was 
developed for the usual, three dimensional Vlasov-Poisson system, i.e., / = 
in the above, in p[10l[IIl[12l[T3l[20], see also Mi [24] . The approach has 
also been used to prove the existence of stable steady states for flat galaxies 
without a halo, i.e., with / = in the above, cf. P, [7J [19]. The fact that 
in the present situation the energy is a functional acting on two functions 
together with the potential interaction terms between the flat and the non- 
flat component requires substantial new ingredients in the basic scheme. One 
pitfall to avoid is that for a minimizer of the above energy functional one of 
the two components might vanish. 

Besides the above stability results it is known that global classical so- 
lutions to the initial value problem for the usual three dimensional Vlasov- 
Poisson system exist, cf. [HI [18], while local classical and global weak solu- 
tions exist in the flat case without halo, cf. [4]. For the situation at hand 
nothing is known about the initial value problem, but we conjecture that 
the analogue of [4] for weak solutions remains true. Our stability result is 
conditional in the sense that it holds for solutions as long as they exist and 
preserve the required conserved quantities. For more information on the 
Vlasov-Poisson system in general we refer to the review article (22]. 

The paper proceeds as follows. In the next section we formulate our 
variational problem and our main result on the existence of minimizers. In 
Section [3] we establish properties of the potentials which allow us to define 
and control the potential energies, in particular the interaction terms. Next 
we collect some relevant results about the decoupled variational problems 
where one of the two components is missing; these facts are established in 
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an appendix. In Section we show that the energy functional is bounded 
from below, that not all the mass can escape to infinity along a minimizing 
sequence, and we investigate the splitting properties of the functional. With 
these prerequisites we can then prove the existence of minimizers in Section [6j 
The fact that such minimizers are steady states together with some of their 
properties are established in Section [7J In Section [H] we finally investigate 
the stability estimate resulting from their minimizing property. 



2 Variational setup 

We denote the set of non-negative, Lebesgue integrable functions by L^_(IR n ). 
For / e L\ (M 6 ) and peL\ (R 3 ) we denote by 

p f (x):= [ f{t,x,v)dv, U p {x):=- [ P ^ dy 

the induced spatial density and gravitational potential; we write Uf = U Pf . 
Similarly, for /6L^(1 4 ) and peL^(M 2 ), 

ps(x):= [ f(t,x,v)dv, Up{x):=- [ ^ f„ dy, 

and to abbreviate we sometimes write p and U instead of pj and Up] notice 
that the latter is defined on IR 3 . In what follows we do not explicitly denote 
the domain of integration — IR 3 or M 2 — unless in cases of ambiguity. The 
integrability properties of these potentials are investigated in Section [31 Next 
we define the various parts of the energy functional. For /6L^(M 6 ) and 

£odn(/):=^ J J \v\ 2 f(x,v)dvdx, E kin (f):=^ J J \v\ 2 f(x,v)dvdx, 

Bpol (/):=-i ffei&M**, Bpot (/ ): =-i ff^piMdxdy, 
2JJ \x-y\ 2JJ \x-y\ 

denote the kinetic and potential energies of the non-flat and flat components. 
The total energy of each component is then defined by 

H(f): = E kin (f)+E pot (f), n(f): = E kin (f)+E pot (f)- 
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Finally, 

n{f,f)=H(f)+H{f) + ± J U f (x)p f (x)dx + ± J U f (x,0)p f (x)dx, 

= H(f)+H{f) + J U f ~(x)p f (x)dx 

is the total energy of the state (/,/). In Section [3l where we investigate the 
existence of all these integrals on the constraint set defined below, we will 
also see that the two interaction terms are equal. 

We wish to minimize this functional over the constraint set 

-^M :={(/,/)! /g£+(IR 6 ), /eL+(M 4 ), ||/||i<M, ||/||i+i/fc<iV, 
I l/l |i < M, 1 1/|| 1+1/ * < N, E kiQ (f) + EUf) < oo, 
f(x,x 3 ,v,v 3 )=f(x,-x 3 ,v,-v 3 )j, 

where M := (M, N,M, N) denotes the constraint vector whose components 
are all strictly positive, 1 1 • | \ p denotes the usual L p norm, and 

0<k<7/2, 0<~k<2. 

In Sections |3] and [5] we will see that the total energy functional is well de- 
fined and bounded from below on this set. The constraints on ||/|| 1+1 / fc and 
I l/l li+i/fc P^y the role of the Casimir constraints, and it does not seem to 
be possible to include these Casimirs into the functional to be minimized, as 
was done for example in [11] for the purely three dimensional and in [19] for 
the purely flat problem. The following theorem is our main result. 

Theorem 2.1. Let (fj,fj) C Tyi be a minimizing sequence oj It. Then 
there exists (/o,/o)^^M) a subsequence again denoted by (fj,fj) and a se- 
quence of shift vectors (o 3 -) CM 2 such that with Tjfj(x,v):—fj(x-\-(dj,0),v) ! 
Tjfjix^-^fjix + a^v), 

Tjfj^fo, Tjfj^fo weakly in L 1+1/fc (M 6 ) orL 1+1/ ^(M 4 ) respectively, 
Epot (Tj fj - f ) — > , E pot (Tj fj - f ) — > , 

and 

J (/'if -Pf.)'', , ///•'• -('• 
Moreover (/o,/o) is a minimizer of It over T^j\_. 
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The spatial shifts parallel to the (xi,x 2 ) plane are necessary due to the 
invariance of the total energy and the constraint set under such shifts. If 
(/o,/o) is a minimizer of H, then (Tjf ,Tjf ) is a minimizing sequence for 
any choice of shift vectors dj G M 2 which is weakly convergent to a minimizer 
only if we shift our frame of reference accordingly. 



3 Preliminaries 

We start by collecting some well known estimates for the spatial densities 
and potential energies induced by elements from the constraint set JF M . 

Lemma 3.1. Let (/,/) GJF M and define n := k + 3 / '2 , h:=k + l. ThenpfE 
L 1+1 / n (R 3 ), p f eL 1+1 ^(R 2 ) with 

\\Pf\\i + i/n < CN^^EUf) 1/Ck+2 \ 

and 

-E pot (f)<C\\Pf\\l/ 5 <CME km (f) 1/2 , 

— Epot(f) < C||P/|l4/3-CM-Ekin(/) 1/2 ) 

where the constant C > is universal and Cm > depends on the constraint 
vector M. By the restrictions on k and k, l + l/n>6/5 and l + l/n>4/3 so 
that Pf eL 6 / 5 (R 3 ), Pf -eL 4 / 3 (R 2 ). 

Proof. Given R > we split the ^-integral and use Holder's inequality and 
the definition of the kinetic energy to find that 



Pf( x ) = / f(x,v)dv+ / f(x,v)dv 

J\v\<R J\v\>R 

(A \ l/(fc+l) / /• \ fc/(*+l) i rr 

< l^—R 3 j U f 1+1 / k (x,v)dvj + — JJ \v\ 2 f(x,v)dvdx. 

We optimize this estimate in R, take the resulting estimate for Pf(x) to the 
power 1 + 1/ n and integrate with respect to x to obtain the estimate for pf. 
The estimate for pj follows the same lines. The last two inequalities follow 
by interpolation and the Hardy-Littlewood-Sobolev inequality. □ 
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In order to analyze the mixed term in TC(f,f) we need some information 
on the integrabihty of the flat potential in R 3 . 

Lemma 3.2. Let pGL 4/3 (M 2 ). Then UpeL 6 (R 3 ) and 

H^pIIl 6 (r3) <C||pIIl 4 /3(ir2). 

Proof. We use the general form of the Minkowski inequality, cf. [T6l 2.4], and 
the weak Young inequality to obtain 



\Up\ 



P(y) 



< 



r2Jr W k-(y,o)| 
p 6 (y) 



dy I dx 3 dx 



C 



{\x-y\ 2 + xlf 

6 



1/6 • 

dx 3 ) dy 



dx 



p(y) 



\x— y\ 5 ^ 6 



dy 



^ = C||5*|-|- 5/6 || 6 r 



^ C||P|IL 4 /3(R2)| 



. 1-5/61 |6 



12 / 5 fTn>2^ 



the function 
4.3]. 



is an element of the weak LP space C /A (M n ), cf. pi 



□ 



We also need to investigate the integrabihty of U p , restricted to the (xi,X2) 
plane. 

Lemma 3.3. Let pGL 6/5 (M 3 ). Then U p (-,0) G L 4 (R 2 ) with 

\\U(',0)\\ L 4<p!>)<C\\p\ \ L 6/ B{w3 y 

If in addition p E L^ 3 (R 2 ) , then the following mixed potential energies exist 
and are equal: 

Up(x)p(x)dx= / U p (x,Q)p(x)dx. 



Proof. Fubini's theorem together with the Holder inequality and Lemma I3T21 
imply that for pGL 6/5 (K 3 ), pGL 4/3 (M 2 ), 



UJx 



r,0)p(x)|dx< f[ \^ y l& ) \ dxdy= [ \U m (y)p(y)\dy 



JJ \(x,0)-y 
<C\\ P \ 



z, 6 / 5 (r 3 )IIPIIl 4 /3(r2). 
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The estimate for ||C/(-,0) ||l4(k 2 ) follows by taking the supremum over all 
pGL 4 / 3 (R 2 ) with | |p| | L 4/3( R 2) = 1. Since the mixed potential energies now 
exist they are equal again by Fubini's theorem. □ 

It will be useful to view the potential energy as a bilinear form which 
induces a scalar product. More precisely we define for 

If p(x)a{y) 
(p,0-) P ot:=7T / -j r dydx 



)3> 



2 J \x — y\ 

with the analogous definition for (p,a) pot , p,a G L^ 3 (R 2 ), and 

It is well known that (v)pot is a scalar product on L 6 / 5 (R 3 ), cf. [HI 9.8], and 
the same is true on L 4 /3(M 2 ). The induced norms are denoted by 

||p|| P ot:=(p,p)pot ) I |p|| P ot := <P, ^pot- 
Finally, ( f,g) pot '■= (Pf,p g )pot etc, provided that the induced spatial densities 
belong to the proper LP space, so that with this notation, 

^ P ot(/) = -(/,/)pot = -||/||^ (3-2) 

etc. The Cauchy-Schwarz inequality corresponding to the mixed case (13.11) 
is established next. It tells us how strong the mixed potential energy term is 
in comparison to the potential energies of its individual components. 

Lemma 3.4. Let p e L 6 / 5 (R 3 ) , peL^ 3 (R 2 ). Then 

|(P,P)pot|<||p||pot||p||pot- 

Proof. We first show the assertion under the additional assumption that 
p, p G are compactly supported and smooth. In that case U p is smooth 
and bounded. Let deC™(R 3 ) be such that d> and fd= 1, and let 5 e (x) : = 
e~ 3 d(x/e) denote the induced 5-sequence; e>0. Then 



lim / U p (x,xs)6 e (x3)dx3 = U p (x,0) 
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pointwise for iel 2 , while the latter integral is bounded in modulus by 
H^plloo- Using Lebesgue's theorem and the fact that (v) po t is a scalar prod- 
uct on L 6 / 5 (IR 3 ) we can now argue as follows: 



UP) 



pot I 



U p (x,0)p(x)dx 



lim|(p,p<g>5% ot | 

e— +0 



< Hpllpotlim - 







= -lim 


// 







U p (x, x 3 ) 5 e (x 3 )p(x ) dx 3 dx 



< 



lim ( - 



pot 



\P \ Ipot I \P \ Ipot) 



p{x)p{y)5 e {x 3 )5 e {y 3 ) 
\x-y\ 

p(x)p(y)5 e (x 3 )5 e (y 3 ) 
\x-y\ 



■dxdy 



dxdy 



1/2 



1/2 



in the last step we used that J 5 e = 1 for e > 0. The general case follows by 
approximating p and p by compactly supported, smooth functions, observing 
the fact that both sides of the inequality are continuous with respect to the 
L 6 / 5 (IR 3 )-norm for p and the L 4 / 3 (R 2 )-norm for p. □ 



4 The decoupled minimizers 

In the next sections the existence and properties of the minimizers of the 
decoupled problems where one of the components is missing will become 
important. Here we briefly collect the relevant facts. A function g on M. d x M. d 
is called spherically symmetric if for every AgSO(gT), g(Ax, Av ) = g(x,v ). 
For each M,N > the energy Ti.(f) has a minimizer /q D in the set 

^ N :={feLl(R 6 )\ ||/||i<M,||/|| 1+1/Jfc <JV, £ ki n(/)<oo}. 

The minimizer is unique up to spatial shifts, spherically symmetric, has neg- 
ative energy, i.e., 7i(/o D )<0, saturates the constraints, i.e., ||/q D ||i = M, 
| |/o D | |i+i/fe = N, and has compact spatial support. There exists a constant 
R* > which is independent of M and N such that the radius of this spatial 
support is 

R = R * M (2k~l)/3 N -(2k+2)/3^ ( 4 _-Q 

By spherical symmetry, 

/o D &x 3 ,v,v 3 ) = / 3D (x,-x 3 ,v, -v 3 ). 
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Similarly, the energy Tt(f) has a minimizer f^ L in the set 

Fm,n '= { f G L\ (M 4 ) | ||/||i<M, \\f\\ 1+1/ - k <N, E kiD (f) < oo } . 

A slight complication arises from the fact that we do at the moment not know 
whether this minimizer is again unique up to spatial shifts. However, there 
does exist a two-parameter family (fJ^ N )M,N>o such that ff^ N is a minimizer 
of Tl(f) over J-J} J N which saturates the constraints, has negative energy, is 
axially symmetric with respect to the £3 axis, i.e., spherically symmetric as 
a function of x,v, and has compact spatial support. There exists a constant 
R* independent of M and N such that the radius of this spatial support is 

R = R*M~ k N- { ~ k+1) . (4.2) 

In what follows / FL always denotes the corresponding member of the above 
family. In particular, if 

M = (M,N,M,N)=: (M 3D ,M FL ) 

then /q D denotes the minimizer of Ji over T^ N and f$ h denotes f^^- 

Since the above facts are known or follow by arguments already available 
in the literature we defer their discussion to the appendix. 

5 Properties of H 

First we establish a lower bound for H on JF M and certain a-priori bounds 
along minimizing sequences. 

Lemma 5.1. (a) The functional Ti is bounded from below on JF M; i.e., 

— 00 < inf H =: h M < 0. 

(b) Along every minimizing sequence (fj,fj) C J~m of Ti both the kinetic 
and the potential energies are bounded, more precisely, for j sufficiently 
large, 

Ekm(fj)+E kin (fj) + \E pot (fj)\ + \E pot (fj)\ <Cm, 
where the constant Cm>0 depends only on M. 
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Proof. Lemma [3.11 and Lemma [3.41 imply that for (/,/) GJF M , 

CfJ>pot| < ll/llpotll/llpot^CM^kinl/) 174 ^^) 174 
<C M ^kin(/) 1/2 + C M ^kin(/) 1/2 . 

Using Lemma 13.11 again this yields the estimate 

H(fJ)>EUf)-CMEUf) 1/2 + EUf)-CMEUf) 1/2 - (5-1) 
Hence /^m > — oo. Moreover, 

ftM<Wo 3D ,/o FL ) = «(/o D )+Wo FL ) + J UoPodx<0. 
Hence along a minimizing sequence T~t(fj,fj) <0 for j sufficiently large, and 

byflEU), 

{E kin U)) 1/2 - Cm/2) 2 + (E kin (fj) 1/2 ~ Cm/2) * < C M /2. 

Another reference to Lemma 13.11 completes the proof. □ 

In order to pass to the limit along a minimizing sequence we need the 
following compactness properties of the potential energies; by Is we denote 
the indicator function of the set S, and we recall ( 13.21) and the corresponding 
notation. 

Lemma 5.2. Let ( Pj ) C L 1+1 / n (M 3 ) and (pj) C L 1+1 /"(M 2 ) be such that 

pj^po weakly in L 1+1/n (R 3 ), pj ->> p Q weakly in L 1+1/h {R 2 ). 

Then for each R>0, 

ll 1 Si i (Pj-Po)||pot-^0, ||lB J? (Pj-po)|| P ot-^0 asj^oo. 

Proof. The convergence of the non-flat potential energy is proved for example 
in [22l Lemma 2.5]. For the flat case we refer to [7J Lemma 3.6]. □ 

A crucial step in the analysis is to show that minimizing sequences do not 
spread out in space and that up to spatial shifts not all the mass can leak 
out to infinity. This is the content of the next result. 
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Proposition 5.3. Let {fj,fj) C J-"m a minimizing sequence of 7i. Then 
there exists a sequence (a,) C M 2 0/ shift vectors, e > 0, and R >0 such that 
for all sufficiently large j e N 7 

/ / fjdvdx>e , / / fjdvdx>e . 

JiajM+BRgJ Jdj+B Ro J 

Here B Ro and B Ro denote the closed ball of radius R about the origin in M 3 
or M? respectively. 

Remark. It is important that the same shift vectors work for both the 
non-flat and the flat component. 

Proof. Let Uj :—Uf v pj '—P^, and let R 3D and R Fh denote the radii of the 
decoupled minimizers /q D and f FL subject to constraints M 3D and M FL , cf. 
Section HI Since 

limH(f j J j )<H(f™J Fh ) = H(f™)+H(f^) + [ U™p F L di, 
we get that for j sufficiently large, 

/ uf Pj dx<n(f^)+H(f FL ) + ^ j u^ D p FL dx. 

Since H(fj) > H(f$ D ) and H(fj) > H(f£ L ) this implies that 

/ Uf Pj dx < \ J Ufpl-d~x = -\ JJ ^Iff^ dxdy 
MM 

< ~2(R™ + R^) (5 ' 2) 

for all sufficiently large j e N. 
For R > 1 we write 

11 11 

TZT = 1 {|x|<i/fl}W]-| +Mi/R<\x\<R}(x) — + l { \ x \> R} ( y x) — 

=:K R (x)+Kl(x)+K R (x). 
With this splitting 

Pj(y)Pj(Z) 



UjPjdx 



\{x,o)-y\ 
12 



dx dy = J\ + J2 + J3 



The second and third terms are estimated straightforwardly: 



J2<R p j (y)p j (x)dxdy, 

JJ\(x,0)-y\<R 

J3 < iT 1 // p j {y)p j {x)dxdy<MMR- 1 . 



For J\ we first apply the Holder inequality and then the general form of the 
Minkowski inequality as in the proof of Lemma 13.21 to obtain the estimate 



J\ < HPjHl+l/r 



Pj{x) 



-dx 



i-~\<l/R 1(^,0) 

<C||P,l|l+l/n||p,*(^) n/(n+1) ||n + l 
— C||Pjl|l+l/n||Pj||l+l/n-R CT <C*M-R 



n+1 



where 



7: 



1 1 

+ - 



>1, er:=-- 



n 



>0; 



^n + 1 n + 1 / 7 n+1 

recall that 3/2 < n < 5 and 1 < n < 3. With fl5\2l we find that 

MM 



and hence 



2(R 3D + R FL ) 

MM 



> J Ujpjdx = — Ji — J2 — J3 



MM 



— CmR a ■ 



2{R 3D + R FL ) R 
For i? sufficiently large the right hand side is positive, so that 



(XiT 1 



MM 



MM 



2{R™ + R FL ) R 



x-(y,0)\<R 



C M R " ] < // p j (x)p j (y)dxdy. 

(5.3) 

The existence of the shift vectors (a,,) with the asserted properties is now a 
consequence of the following lemma. □ 

Lemma 5.4. Let peL\(R 3 ), aeL\(M. 2 ) with 

0< / p(x)dx, / a(y)dy<M < 00 
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and such that for some So, r > 

// p(x)a(y)dxdy> Sq. 

J J\x-(y,0)\<r o 

Then there exist 6q, Rq>0 depending only on Sq, ro, and M such that 
e < / p(x)dx and e < / a(y)dy 

J\x-(a,0)\<R J\y-a\<R 

for some a£l 2 . 



Proof. Let zeM be given. Note first that 

{(x,y)GM 5 | \z-(y,0)\<r , \x - (y,0)\ <r } 

C {(x,y)eM 5 | \z-x\ < 2r , \x - {y, 0) | < r } 

and hence 

a{y)dx J dy 



Z-{y,0)\<r V \x-(y,0)\<r 



< / I / cr(y)dy\dx. 

J\z~x\<2r \J\x-(y,0)\<r J 

Multiplying with p(z) and integrating with respect to zel 3 we obtain 
p(z) / I / a(y)dx J dydz 

J\z-(y,0)\<r \J\x-(y,0)\<r J 



<[p(z)[ (f 

J J\z-x\<2r \J\x- 



o~(y)dy ) dxdz. (5.4) 

<y,o)\<r 



Changing the order of integration, the right hand side of (j5.4p can be rewrit- 
ten as 



/ / p{z) ( / <r(y)dy) dzdx 

Jx^M. 3 J\z-x\<2r \J\x-(y,0)\<r / 



p(z)dz ] I / cr(y)dy J dx 

x&R 3 \J\z-x\<2r J \J\x-(y,0)\<r 

R(x)S(x)dx, 
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where the functions R and S are defined by 



R(x):= p(z)dz, S(x):= v(y)dy. 

'\x-z\<2r J\x-(y,0)\<r 



From our hypothesis, 



5 < I p(z) I a(y)dydz 

\z—{y,0)\<ro 



T— 3 / P[z) [ [ cr(y)dxdydz, 

47rr JzGM. 3 J\z-(v,0)\<rn J\x-(vM<m 



JzGM. 3 J \z-(yfi)\<r J \x-(y,0)\<r 

so that combining with (15.41) we are led to 

5 < q / R(x)S(x)dx. 

As a direct consequence of our definitions (I5.5P we find that 

\\R\\oo<M, ||^||oo<M. 

Furthermore, 



= / / p(z)dzdx= / / p(z)dxdz 

JxGR 3 J\z-x\<2r JzeR 3 J \z-x\<2r 

4:71 

< 8M— r 3 , 



and 



\\S\\i= / a(y)dydx= / a(y)dxdy 

JxeR 3 J \x-(y,0)\<ro J yeM. 2 J \x-(y,0)\<r 

, Ait o 

3 

We may thus continue with (15. 6p as follows: 

3 

"o 



5q < o / R(x)S(x)dx 



< -l-MRS) 1 ' 2 ^ I {R{x)S{x)) l ' 2 dx 



4tiTq 



i/ 2 / r \ 1 / 2 

-llff^ll 1 / 2 ( / TK^dT. \ 
< 2v / 2M|| J RS||^ 2 . 



<l^\\ RS \\)L 2 (f R{*)dx) [ I S(x)dx 



15 



So there exists a G M 3 such that 

R(a) S(a) > 
In view of (15.71) this implies that 



<5n 



V 



2 V / 2M/ 



i*(a) > and S(a) > ^ 



]M 3 v ' 8M 3 ' 

Finally we write a = (0,03) with a G M 2 , 03 G K and observe that 



/ cr(y)dy> / a (y) dy = S (a) > 

J\a-y\<ro J \(a,a.3)-(y,0)\<ro 

In addition, S(a) >0 clearly implies |a 3 | <r , so that 

/ P(z)dz> [ p(z)dz = R(a)>-^ 

J \z-{a,0)\<3r J\a-z\<2r ° 1V1 

which is exactly our claim with 



2 


M 3 ' 



e ° := 8M3-' ^o:=3r . (5.8) 

□ 

In what follows it is important to control the parameters eo and Rq in 
Proposition 15.31 if the constraint vector M varies. This is the content of the 
following corollary. 

Corollary 5.5. Let the constraint vector M satisfy the bounds 
()-. .\/,< M < M„. 0<Mi<M <M UJ 

0<Nt<N<N u , 0<Ni<N<N u . 

Then the parameters eo and Rq in Proposition 15.31 can be chosen indepen- 
dently o/M and (fj,fj), depending only on the bounds Mj and M tt . 

Proof. Under the given bounds on M we can choose R > depending only on 
these bounds such that the left hand side in (15.31) is bounded from below by a 
parameter <5 > also depending only on these bounds. To this end, observe 
that M and M are bounded both from below and above, Cm is bounded 
from above, and the radii R 3T> and R FL are bounded from above in view of 
P~T|) and Given (g£gj) this completes the proof. □ 
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The last tool needed for the proof of Theorem 12.11 is the fact that the 
energy infimum h-^ is sub-additive in M. While up to now all components 
of the constraint vector M were strictly positive, this sub-additivity is for 
technical reasons needed also in situations where the flat or the non-flat 
component of a constraint vector vanishes, i.e., M= (M 3D ,M FL ) and M FL = 
or M 3D = 0. In such a case /im is obviously taken to denote ?i(/o D ) or 
7i(/ FL ) respectively, where /q D is the minimizer of H over J-'^u and / FL 
is the one over ^Jpl , cf . Section HI We say that the constraint vector M G 
[0,oo [ 4 is nontrivial, if 

(M>0AiV>0)V(M>0AiV>0). 

Proposition 5.6. For all Mi, M 2 G [0,oo[ 4 ; 

If both Mi andM 2 are nontrivial, then this inequality is strict. //Mi satisfies 
uniform bounds from above and below as in Corollary \5.5\ and if either this 
is true also for M 2 or one component of M 2 vanishes and the other one 
satisfies such uniform bounds, then there exists e > depending only on these 
bounds such that 

^Mi+M 2 < ^Mi + h-M 2 — e- 

Proof. Consider two minimizing sequences (fj,fj) CJF Ml and (fj,ff) C Jm 2 
with 

7~i{fj Jj)^h Ml , T~L(ff,ff)^h M2 . 

If one of the constraint vectors, say M 2 has a trivial component, say the flat 
one, then we take (/q D ,0) as the corresponding minimizing sequence. By the 
Minkowski inequality, (fj + ff,fj + ff)EJ r M 1 +M 2 , an d hence 

^Mi+M 2 < Ti{f] + fjjj + fj) 

Hifj.f!) ■H(fj.fj) 

~ ) fj )pot ^ifj 1 fj )pot 2 {fj , fj ) po t 2(/y , fj )p t 

<H(f},f})+H(ff,ff)^h Ml +h M2 . 

If Mi and M 2 have both at least one nontrivial component, then the corre- 
sponding potential interaction energy is strictly negative so that the estimate 
above is strict. 
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Assume now that we have positive uniform lower and upper bounds for 
Mi and M 2 . We can assume that the minimizing sequences are shifted in 
such a way that the assertions of Proposition 15.31 hold with ej, 6q, R^, and Rq, 
and without spatial shifts. If one component of M2 vanishes, say, the flat 
one, the corresponding trivial minimizing sequence (/o D ,0) need of course 
not be shifted, and we take for R% the radius of the minimizer /q D and for 
6n its mass. In either case 



1^2 



(p>f) P ot> // —. ^-r-dxdy^-^ 



o c o 



"o 

The latter quantity is now bounded from below by some e > depending only 
on the uniform bounds on the constraint vectors, where we use Corollary 15.51 
and in addition (14. ip if the flat component of M2 vanishes. If the non-flat 
component of M2 vanishes we use (14.21) instead. □ 

Remark. The uniform sub-additivity is also valid if both Mi and M 2 have 
exactly one nontrivial component which is uniformly bounded from below 
and above, but this case is not needed in what follows. 



6 Proof of Theorem 2.1 



Let (fj,fj) G J-~m be a minimizing sequence for Ti. We choose shift vectors 
fij G M 2 such that the assertion of Proposition 15.31 holds. To keep the notation 
simple we redefine {fj,fj) as the minimizing sequence shifted by these vectors 
as in the statement of Theorem 12.11 Hence according to Proposition 15.31 

eo< / fjdvdx<M, e < / fjdddx<M. (6.1) 

J\x\<R J J\x\<R J 

This new sequence is of course minimizing as well. The definition of JF M 
implies the a-priori bounds 

||/illi + i/ fe <iV, \\f 3 \\ 1+irk <N. 

Hence after extracting a subsequence which we denote by the same symbol, 

f 3 - f weakly in L 1+1 / fe (M 6 ), f 3 - f weakly in L 1+1 /\r 4 ). 
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From this weak convergence it follows that 

||/o||i<M, ||/o||i<M, ||/o||i + i/ fc <iV, ||/o|| 1+1/ fc<iV, 

and 

Ekm(fo) < limsupE kin (/ i ) < oo, E kin (f ) < limsup£ kin (/j) < oo. 

By Lemma [XT] the corresponding spatial densities Pj\—pf j and pj: = p^ are 
bounded in L 1+1 / n (IR 3 ) or L 1+1 ^ n (E?) respectively. After extracting a subse- 
quence again, 

pj -± p weakly in L 1+1/n (M 3 ), pj -± p weakly in L 1+1/fl (R 2 ). 

It is easy to see that in fact Po = P/ an d Po = P/ - The essential step is to 
prove that up to extracting yet another subsequence the potential energy 
terms converge, i.e., 

I \fj - /ollpot + \\fj~ fo\ | P ot -> as j -> oo; 

by Lemma 13^0 it then follows that also (fj — fo,fj — /o) po t - >■ 0. 

For R>R 1 >R we define B Rl>R := {x G M 3 | i?i < |x| <i?} with the obvi- 
ous definition of B Rl R , and we split the functions fj and as follows: 

fj = 1 B Rl xR^fj + 'i-B RltR xRsfj + 'i-B Ri00 xRsf j =: fj + ff + ff, 
fj = 1 B J?1 xR 2 /j + iB^^xK 2 /?' + ^b^xR 2 /? =: fj + fj + fj- 

Lemma 15.21 implies that for R > Ri > R fixed, 

li/i+/|-/o-/ollpot+||^ 1 + ^-/o 1 -/ollpot^Oasi^oo. (6.2) 
So we only need to show that for any e > and R sufficiently large, 

Uminf f||/f|| pot + ||^|| pat )<e. (6.3) 

Once this is established we use the triangle inequality for || • || pot to conclude 
that 

Wfj — fo 1 1 pot < | \fj + fj — fo — fo I Ipot + \\fj\ |pot + \\fo ||pot- 
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We can surely find R > 1 such that the right hand side is as small as we want 
for j sufficiently large. Hence for j — > oo, 

E po t(fj) ^E pot (fo) 

and with the same argument, 

E P ot(fj) — >E pot (f ). 

Finally by Lemma 13.41 

UjPjdx-^ J UoPodx, 

and all together implies that 

H(fo,fa)<limH(f j J j ) = h m . 

This is the desired minimizing property of (/o,/o)- 

We prove (16.31) by contradiction, so assume that (16. 3p is false, i.e. 

3et > OWR > 1 3j e NVj >j : | \ff\ | pot + 1 1^| | pot > 6l . 

Then we can choose a subsequence such that without change of labeling it 
satisfies either 

Vi?>13j GNVj>j : [I//I|pot>ei/2 (6.4) 

or 

VR>13j eWj>j : ||/ J 3 ||pot>ei/2. (6.5) 

In the following we consider the first case, the second one can be treated 
analogously. The contradiction is arrived at by splitting fj and fj as above 
and then using the uniform sub-additivity from Proposition 15.61 Let us 
denote 

fj '■= ^B Ro xR 3 fj, fj '■= l^xR 2 /?'- 

Since the splitting parameters satisfy the relation R> Ri>R , (16. ip implies 
that 

eo<i|/ J ||i<ii//||i<M, e <||£%<ii//|ii<M. (6.6) 

We also need uniform lower bounds for the L 1+1//fc -norm and L 1+1 / fc -norm. 
By Lemma [3.11 

\\f% = \\P% < C(Ro) \\ P % + yn < C(Ro) \\f°\t + $ (n+1 \ 
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with an analogous estimate for fj. Hence with (16. ip . 

< C(e ) < \\fj\\i + x/k <N,0< C(e ) < \\fj\\ 1+lfk < N. (6.7) 

From the assumption (16. 4p we now derive such bounds also for By 
Lemma 13.11 with #g]0,1[ an interpolation parameter and er:=(l — #)(1 + 
fc)/(l + n), 

II^IIS^ = <C|| / ,||i /5 ^CII^II^II^H^^ <<C||^||^||y-||?- ;L/fc . (6.8) 

With / = fj this implies that 

0<C(ei)<||/?|| a <M, 0<C{e 1 )<\\ff\\ 1+1/k <N. (6.9) 

To arrive at a contradiction we insert the splitting of fj and fj into the energy 
functional: 

mfrf,) = mfj.fl) ■H(fj.fj) ■mjj.r;) 

~2(fj)fj+fj)pot — 2(fj,fj)pot — 2(fj,fj+fj)pot — 2(fj,fj)pot 

~^(fjlfj + fj )pot — 2(/j , fj) pot ~2(fj + fjl fj) pot ~~ fj) po t 

^mfj.f*) -mfj.fj) -mjj.f;) 

-h -I 2 - h -Iz-J.-Jz-J,- J 2 . (6.10) 

Using the Cauchy-Schwarz inequality for (v)pot, i-e., Lemma I3~4l for the 
mixed terms Ji and Ji, and the boundedness of potential energies along the 
minimizing sequence, cf. Lemma I5.1[ we obtain the estimates 

h + Ji<C\\ff\\ pot <C{\\ff-fXot + \\fo\\pot), 

h + Ji < c i | pot < c (| \ff - ; 2 | u + 1 l/ 2 l | P ot) . 

For R > 2R 1 and x G -Bj? 1 , y G .Br,oo we note that 

1 1 1 _ 2 

and we combine this with the Holder inequality to estimate I2, h, J2, and 
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J 2 as follows: 

1/2 



l2 - 2 f B Pj(%)dxJ^ \y\ 1 Pj{y)dy<C\\p j \\l / s 
/ 2 <2^ pj(x)dxj_ |yrV J (y)^<C'l|Pillv3^^ 



1/2 



J 2 < 2 



/ pj(x)dx f \y\ 1 p j (y)dy<C\\p j \\ e/5 \\p j \\ 4/3 (^jf) 
Jb Rx Jb Ri00 \ n J 



J2<2 pj(x)dx \y\ 1 Pj(y)dy<C\\p j \\ 6 / 5 \\p j \\ 4 / 3 



R 



We wish to apply the uniform sub-additivity from Proposition 15.61 to the 
constraint vectors induced by (//,//) and (ff,fj)- To this end, let 

Mi:=(ll/illi,ll/il| 1+ i/fc,ll/illi,ll/illi +1 /s), 

M} := i=l,2,3. 

From (16.61) . (16.71) . (16.91) we have the required uniform bounds for Mj and 
(Mj,Nj). With respect to ff we now distinguish two cases. Either this 
function also satisfies such non-zero uniform bounds or it is negligible, more 
precisely: 

Casel: 3t 2 > OVi? > 13j E NVj > j : \E pot {ff)\ >e 2 . 
In this case the analogue of the potential energy estimate (16. 8p for /? implies 
that 

< C(e 2 ) < 1 1/7| U < M, < C(e 2 ) < \\ff\\ 1+l/ ~ k < N. 

So we have obtained uniform positive bounds for each entry of the quantities 
Mj and M| from above and below. By Proposition 15. 6[ 

n(fjJI)+n(ffjf)+n(ffJf)>h M i+h m 2+h M s>h M +e 

with e > independent of the splitting parameters R > 2Ri and of j. Recall- 
ing (16.101) we find that 



hM-n(f j J j )+e<I 1 + I 2 + I 1 + I 2 + J 1 + J 2 + J 1 + J 2 



I \fo I | pot + I \fo 1 1 pot + I \fj ~ /o I I pot +\\fj~ /o I I pot + {R\/ R) 



1/2 
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We choose R\ > R such that 

Ci(||/o 2 ||pot + ||/o 2 ||pot)<6/4. 
Next we choose R>2R 1 such that d(R/ R^ 2 <e/4. For j large, 



1 



h^-H{f j ,f j ) + e<^e + C 1 \\f'-mU + \\fhfi l 



j JO I Ipot T | \ Jj Jq \ |pot 

and by (16.21) this contradicts the fact that (fj,fj) is minimizing. 

Case 2: Ve > 03R € > 1Vj e N3j >j : \\ff\\ pot < e, provided > # e . 
In this case we neglect /? in the sub-additivity argument and recall that 
Proposition 15.61 yields e 2 > only depending on the bounds for Mj, Mj, and 
iVj such that 

^Mj + + h(Mf,Nf, 0,0) — ^Mj+(M j 3 ,Ar3,0,0) + e 2 + 

> ^M)+(M?,Ar?,0,0)+M? + e 2 

> ^Mj+(M?,iV?,0,0)+M| +^2 + ^(0,0,M|,JV?) — ^M J +C2- 

By the assumption of the present case we can choose a subsequence which we 
keep on denoting as before such that | \f^\ | pot < e for all j G N, where e will be 
determined in terms of e 2 below; if necessary we increase R so that R>R £ . 
By Lemma [3.41 

|<j*/!>pot|<C6. 

Moreover, we choose R>2Ri>2R such that in (16.101) . IH h«/2<e- 

Hence by flBTTOj) . 

>H(f},f}) + H(f]j~])+H(f])-e 2 -Ce 
> h M i + h M 2 + h( M 3 tN s tQtti) - e 2 - Ce 



> h Mj +e 2 -e -Ce. 

If e is chosen properly in terms of e 2 , 

T-t{fj,fj)>h M + e 2 /2 as j^oo. 

This contradicts the minimizing property of (fj,fj)- If one considers the case 
(16.51) instead of (16.41) . all the arguments remain the same with the roles of 
the flat and non-flat components interchanged. The proof of Theorem 12.11 is 
complete. 
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7 Properties of the minimizer 

First we exclude the possibility that for a minimizer / = or / = 0. Indeed 
the next result shows that the constraints are to some extent saturated by 
any minimizer. 

Proposition 7.1. Let (/o,/o) £ J~m. be a minimizer ofTC over JF M . Then 

||/o||i = MV||/o||i = M, 

||/o||i+i/* = JV, ll/o||i + i/* = JV. 

Proof. We define for a,b,c,d,e>0 a rescaled state (/ *?/o) as 

fo(x,v) :=af (bx,cu), fS(x,v) := df (bx,ev); 

because of the mixed potential energy term x and x must be scaled in the 
same way. Then 

£kin(/o) = ab- 3 c- 5 E kin (f ), E kin (f*) = db- 2 e^E kin (f ), 

E pot (fZ) = a 2 b~ 5 c- 6 E pot (fo), E pot (f*) = d 2 b- 3 e- 4 E pot (f ), 

J UQpldx = adb~ 4 c~ 3 e~ 2 J £7oPo dx. 

Assume that ||/o||i+i/fe < N. Then we choose 

a = c 3 , b = d = e = 1. 

For this choice of parameters f^ — fo, 

ll/olli = ll/olli, H/olli+iA = c3/(fc+1) ll/o||i+i/fe ) 

and 

Wo* Jo) = c~ 2 ^kin(/o)+^pot(/o)+^(/o) + j UoPodx. 

We can choose c > 1 so that the rescaled state still lies in JF M and has 
lower energy which is a contradiction. The analogous argument shows that 
ll/olli + i/i = #- 
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In order to prove that at least one of the two mass constraints is sat- 
urated we assume that ||/ ||i<MA||/o||i<M, and we choose the scaling 
parameters 

a = c -7 , b = d = c~ 4 , e — c. 

For this choice, 

ll/olli = c 2 ||/o||i, ll/olli = c 2 ||/o||i, 
||/olli + iA = c (2fe - 7)/(fc+1) ||/o||i + i/ fc , li/^ll^/^c^V^ii/oli^, 

and 

WoJo)=Wo,/o). 

Since < k <7/2 and 0<&;<2 we can choose c>l such that (fo,fo)£J~M 
and both M^=(||/*|| 1 ,||/*|| 1+1/fc ,||/*|| 1 ,||/*|| 1+1/ ^) and M-M* are non- 
trivial. The strict sub-additivity in Proposition 15.61 implies the desired con- 
tradiction: 

< ^M* + ^M-M* < W(/o )/o) = ^(/o)/o) — ^M- 

□ 

The main result of this section is the fact that the minimizers are func- 
tions of the particle or local energy. We use the Lagrange multiplier method 
presented for example in [TT1 [T2| 122 1 123] . 

Theorem 7.2. Let (/o,/o) be a minimizer as obtained in Theorem \2. II with 
induced potentials (Uo,Uq). Then 

f E -E(x,v) \ k 
fo{x,v)=\^ J a.e., 

f (x,v)=\ ~ I a.e., 

where E(x,v) := h\v \ 2 + U (x) + U (x) and (•)+ denotes the positive part. The 
Lagrange multipliers are defined as 

1 f2k + 5 f \ 

E := j^- (^^— ^kin(/o) + 2^ po t(/o) + J Uopodxj , 

E := -I— (Ck + 2)E kin (fo)+2E pot {f )+ [ U p dx) , 
ll/olli V J J 
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and 

A _ 2(k + l)E kin (fo) - = (k + l)E kin (f ) 

3ii/oii;:;g ' ' ii/oii;;;;* 

Proof. Let (/o,/o) be a minimizer of 7i with corresponding potentials (?7 ;^o)- 
For / such that (/,/o) G.F M we define 

Q(f)-=n(fJo). 

Then 

Q{f)-Q(fo) = E Un {f)-E ]dn (f )+E pat (f)-E pat (f ) 

+ J (pf-po)Uodx. (7.1) 

For each fixed e > we define the set 

Let 77GL°°(R 6 ) be a real- valued function with compact support such that 
?7>0a.e. for (x,v) GM 6 \supp/ and supp^C (IR 6 \supp/ )US' e . FortG[0,T] 
andT = (||7 ? || 1 + ||r ? || 1+1/fe + ||r ? || 0O )- 1 e /2 we define 

/t(x,t;)==a 3 (01|/ol|i ,|^| Mt)v), 
ll/o + ^||i 

where 

m _( \\M\i + i/k n/o+^iii y* +i)/3 

a{> ' V ll/olli ll/o+N|i + i/J 

For tG[0,T], 

||/*||i = ||/o||i, ||/*||i + i/* = ||/o||i+i/ fc 
and fo + tr)>0 a.e. For e small enough, 

n||l <||/o + ^||i<||/o||i + J, 



2 - |UU " |1 - | " uni 2 

^^^ll/o + ^llH-V^II/ollw + l, 
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which implies that a is a smooth function on [0,T] and 



a/(t) = *±I a (i) 



II (fo + tv) 1/k vdxdv 



ii/o+^n; +iA 



n+i/fc 



\\fo+tv\\i 

Moreover, a" is bounded on [0,T]. From (17.11) we conclude that for t G [0,T], 

'"IU , \ 77! / £ \ i Il/Olll* 



Q{f t )-Q{h) 



« 2 (t)||/o + ^lh 



"I 4n(/o) + 



« 2 (*)||/o + t»7ll 



-E kin (r]) 



oil? 



ll/o + ^|| 2 

2+2 



i'- 



+ 
+ 

+ ll/o + ^||i 
By Taylor expansion at i = 0, 

L = -t 



r r pot(/o)+ E/^M/ 



ll/o + NI? 



\\fo+tv\\i 



ll/o + ^lh 



2 ; Prpodx 
poUadx 



p v Uodx. 



oi 



^Wll/o+^lli 

/ 01 II* 



jfc + 1 



oili 



0111 



i IS fi /k ydxdv \ 
o\\i+i/k /J 



(7.2) 



0(t 2 ), 



« 2 (*)ll/o + ^|h 
oil? . 



*+o(o, 

2|Nli* 



ll/o + ^H 2 



■0(t 2 ), 



Olll* 



o M i 

2\ 



ll/o + tr/H 2 
l/o||i i 



t+o(0, 



l/o+tylh 

ll/olli* 



r + 0(* 2 ), 



Olll 

2\ 



t + 0(t z ). 



\\fo+tv\\i 

If we substitute these expansions into ( 17.21) . we find that 

GUt)-G(h)=t fj (E-E + \ti /k ) V dvdx + O(t 2 ) 

with E'q and A as given in the theorem. Since G(ft) attains its minimum at 
t — 0, the choice of 77 and e — * imply that E — E > on M 6 \ supp /o and 

/ p E\ 

fo = I — I a.e. on supp/ - 



A 
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If we repeat this argument with the roles of flat and non-flat states exchanged, 
i.e., for Q{f) '■ = 7~t(fo 1 f) 1 we obtain the assertion for f . □ 

The previous theorem states that for a minimizer (f Q , fo) both components 
are functions of the local or particle energy in the induced potential {7 ,e = 
Uq + Uq. Since the latter is time-independent, the particle energy is conserved 
along particle orbits, i.e., along the characteristics of the Vlasov equations 
( 11. ID and (11.21) respectively. Hence fo and fo satisfy these equations at least 
formally, and we are justified to refer to (fo,fo) as a steady state of the 
system f ll.ll) - fll.4l) . We do not discuss the regularity of this steady state 
further. However, to conclude this section we want to address the question 
whether these states have spatially compact support. 

Proposition 7.3. Let (fo,fo) be a minimizer as obtained in Theorem \2.1\ 
and assume that 

0<k<5/2 and 0<k<l. 
ThenUo, U Q , p Q , PqGL 00 ^ 3 ) with 



lim Uq(x) = 0, lim U (x) = 0, 



E ,E <0, an d Po an d Po have compact support. 

Proof. Consider a density p G HL P (IR 2 ). Then t/p e L°°(M 3 ) with 
lim| x |^ 00 [/p(x) = 0, provided p>2. If p is defined on M 3 then the same is 
true provided p>3/2. We prove this assertion for the flat case. Here 

w> . d y+[ . m _ rfy - 



x-(y,0)\<R \ x ~ (2/) 0)1 J\x-(y,0)\>R \ x ~ (lJ^)\ 



l\&-jj\<y/w^ \x-{y,0)\ 



R 



<(7i2(P-2)/(P-l)||n|| , 

ll p H^({l«-5l<V^f}) R ' 

Since this holds for any R > the assertion follows; notice that for R>0 
fixed the first term goes to zero for \x\ — >oo. By the weak Young inequality 
and Lemma U eL 6 (R 3 ) and U (-,0) G L A (R 2 ), and again by the weak 



28 



Young inequality and Lemma E21 U G L 4 (IR 2 ) and U G L 6 (R 3 ). Hence U 0je G 
L 4 (IR 2 ) flL 6 (R 3 ). If we integrate the relations between / , /o, and ?7 ,e from 
Theorem 17.21 with respect to v or v respectively we obtain the relations 

p=c(Eo-U , e )l, p=c(Eo-U , e (.,0))% (7.3) 

where c and c depend on A and k or A and fc respectively. From the inte- 
grability assertions for the potential we conclude that the spatial densities 
have the required integrability provided 6/n>3/2, i.e., n<4 which means 
k < 5/2, and 4/n > 2 i.e., n < 2 which means k < 1. 

It therefore remains to show that i?o < and £o < as claimed; the as- 
sertion on the support of the densities then follows. Assume that E > 
0. Then for \x\ large, p (x) > c(E Q / 2) n which contradicts its integrability, 
and the same argument works for p$. Now assume that Eq = Q. Then 
Po(x) = c(— U 0:e (x)) n , and this again contradicts the integrability of p since 
— Uo, e >C/\x\ for large |x| and C > 0. We prove this for U , the argument for 
Uo being completely analogous. We choose R>0 such that 

po(y)dy=:m>0. 

A<R 

Next we observe that for \y\ <R and \x\ >2R, 

R 



1 



\x-(y,0)\ \x\ 



< 



(\x\-Ry 



If we restrict the convolution integral defining U to the set {\y \ <R} and 
expand the kernel as indicated the assertion on Uq follows. The same argu- 
ment works for Uq so that indeed —Uo^ e >C/\x\ as claimed. If E = then 
Po(x) = c(— Uo :e (x,0)) n which contradicts the integrability of p. Notice that 
under the present assumptions on k and k it follows that n< 3 and n < 2. □ 



8 Stability 

In this section we show how the minimizing property of a minimizer (/o,/o) G 
JF M leads to a stability estimate. Given a second state (/,/) G JF M and de- 
noting the effective potential of the minimizer by Uo, e a simple computation 
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shows that 

n(fJ) = H{f J ) + jj (^\v\ 2 + U , e (x)^ (f-f )(x,v)dvdx 

+ Jj "Q|^| 2 + C/ 0) e(5,0)^ {f-f ){x,v)dvdx 

~ ll/ _ /o||pot _ ll/ _ /o||pot _2 (/ _ /o ) /-/o)pot- 

With E(x,v) = l\v\ 2 + U , e (x) and 

<(/,/), (/o,/o)) := // E(x,v)(f-f )(x,v)dvdx 



+ yy £?(5,0,iJ,0)(/-/o)(5,i;)tZi;ds (8.1) 
we can rewrite this expansion as 

n(fj) = H(f Q j )+d((fj) J (foj )) 

~ 11/ — /oil pot — ll/ _ /o||pot - 2(/-/o,/-/o) P ot- (8.2) 

We need to show that d((f,f),(fo,f ))>0 with equality only if (/,/) = 
(/(bio)- To this end we restrict ourselves to states (/,/) G JF M such that 

l+l/fc_ / A+l/k [iff f fl+l/k_ f fW/k 



/= y /o, y f i+i/k = j / i+i/K , Jf=Jf°>J f + 1 - J /c . ( 8 - 3 ) 

Remark. From a physics point of view a galaxy and its halo are typically 
perturbed by the gravitational field of some distant exterior object. In par- 
ticular, the perturbation will result in a measure preserving redistribution of 
the particles in phase space, and will hence preserve the constraints in ( 18. 3ft . 
On the other hand, the fact that the perturbations lie in JF M means that the 
stars are only shifted within the galactic plane and not perpendicularly to 
it. This is certainly an unphysical restriction. To remove it is a non-trivial 
problem for future research. 

Using (18. 3p and the strict convexity of the function [0, oo [3 ( i— > ( p for p > 1 
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we find that 
d((f,f),(fo,fo)) 



(E-E )(f-f ) + 



A 



+ (E-E )(f-f ) + 



1 + 1/k 
A 



1 + 1/k 



(f 1+1/k - fo +1/k ) 



(f 1+1/ ' k ~ fo +1/k ) 



> 



+ 



(E-E ) + Xf 1 /k ] (f-f ) 

(E-E Q ) + \fi rk ] (/-/o)>0; 



Theorem 17.21 implies that the last expressions are non- negative, and the strict 
convexity implies that equality holds only if (/,/) = (/o,/o)- 

In order to establish a stability result we now wish to apply the above es- 
timates to the time evolution (f(t),f(t)) of a perturbation of (/o,/o)- Clearly, 
we need to require that (/(0),/(0)) G JF M satisfies the constraints (18.31) . More 
importantly, in view of the fact that nothing is known on the initial value 
problem for the system (ll.ll) - (11.4j) . we have to assume that this system has so- 
lutions t\—*(f(t),f(t)) which preserve the total energy, the constraints ( 18. 3ft . 
and (f(t),f(t)) To keep the rest of the discussion simple we further- 

more assume that the minimizer (/o,/o) is unique in JF M up to spatial shifts. 
If the minimizer is up to spatial shifts only isolated with respect to the dis- 
tance measurement used in the stability estimate below, the result remains 
unchanged. If the minimizers are not even isolated one can prove the sta- 
bility of the whole set of minimizers; we refer to [21] for the corresponding 
modifications of the arguments. 

Stability estimate. Assume the minimizer (/o,/o) is unique in JF M . Then 
for every e>0 there exists 5>0 such that for any solution 1 1— > (/ (t) , / '(£)) of 
M.1\) - (T7J\ ) satisfying the above assumptions the following is true: If 

<(/(0),/(0)),(/o,/o)) + ||/(0)-/o||pot+||/(0)-/o||pot<« 

then 

<*((/(*), /(*)) , (/o, /o)) + 1 ~ /ol Ipot + 1 !/(*) - /o| Ipot < e 

up to spatial shifts parallel to the (xi,x 2 ) plane and as long as the solution 
exists. 

We do not call this a theorem because it is not clear that sufficiently regular 
solutions to the initial value problem do indeed exist. Assuming the latter 
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the proof is by contradiction. If the assertion were false, there exists e > 0, 
a sequence of solutions (ti— > (fj(t),fj(t))) and a sequence of times (tj) such 
that for all j G N, 

^((/,(0)J,(0)) ) (/oJo)) + li/i(0)-/o|ipot + ii/,-(0)-/o|ipot<l/j, (8.4) 

but 

^((/i(^),/i(^)),(/o,/o)) + ||/i(ti)-/o||pot + ||/,-(^)-/o||pot >e (8.5) 

regardless of how we shift (fj(tj),fj(tj)) in space. Now (j8.4p and the fact 
that d is non-negative imply that all three terms in (18.41) converge to zero. 
By Lemma I3T41 this is true also for the mixed term (fj(0) — fo,fj(0) — Jo)pot> 
and since Ti is preserved, (18.21) implies that 

H(f j (t j )J j {t j ))=H{f j (0)J j (0))^H(foJa). 

Since (fj(tj),fj(tj)) G^m this means that (fj(tj),fj(tj)) is a minimizing se- 
quence for 7i in JF M . By Theorem 12.11 there exists a subsequence such that 
up to spatial shifts 

Wfjitj) — /o 1 1 pot + \\fj(tj) — /o|| pot ^0; 

at this point the uniqueness assumption for (/o,/o) enters. Again by 
Lemma I3T41 this is true also for the mixed term (fj(tj) — fa,fj(tj) — fo) po t, and 
(ED with (/,/) = implies that also J;(*i)),(/o Jo)) - 

0. Hence all three terms in (18.51) converge to zero, which is a contradiction. 

9 Appendix: Facts on the decoupled mini- 
mizers 

Here we establish the claims on the decoupled variational problems referred 
to in Section HI Several of these claims, in particular for the non-flat case, 
can be found in the literature. 

Existence. 

In most of the previous investigations the existence of minimizers in the 
decoupled cases was not done exactly for the problems stated in Section HI 
Either the Casimir functional, which in the case at hand corresponds to 
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the norm, was included into the functional to be minimized, i.e., an 

energy-Casimir functional instead of the energy was minimized, and only the 
constraint on the L l norm was posed, or the energy was minimized under the 
constraint that the sum of the mass and the Casimir functional was fixed. 
In the three dimensional case the problem with two constraints in the form 
stated in Section H] has been dealt with in [21 El] • We briefly show how the 
method used in [2] can be adapted to the flat case. 

With the help of the Riesz rearrangement inequality [T6J 3.7] and the 
fact that the kinetic energy as well as the constraints are invariant under 
spherically symmetric rearrangements the problem is reduced to minimizing 
Tl(f) where the functions / in the constraint set have the form 

f(x,v)=(p(\x\,\v\), 

with ip: [0,oo[ 2 ^ [0,oo[ non-increasing in each argument. This monotonicity 
implies that for l<q<l + l/k, 

f q (x,v)\x\ 2 \v\ 2 < C t ' f ' <p q (r,s)rsdsdr<C\\f\\ q q , 
Jo Jo 

f(x,d)\i\ 2 \v\ 4 < C [ [ ip(r,s)rs 3 dsdr<CE km (f). 
Jo Jo 

Hence 

JK ' ' ' \\x\ z \v\ 4 , for \v\ >V\\x\), 

where V(|x|)>0 is an arbitrary function and the constant C depends on 
E kin (f), H/Hi, and ll/H 

quantities which are bounded along minimizing 
sequences. The function g induces the spatial density 

rV(\x\) poo 

p g {x) = C\xl 2/q / \v\ l ~ 2/q d\v\ + C\xl 2 / \v\~ 3 d\d\ 

Jo Jv(\x\) 

= C\i\- 2/q V 2 - 2/q (\x\) + C\x\- 2 V 2 (\x\). 

The choice 

V(\x\) = V q (\x\):=\x\ {1 - q)/{2q - 1) 

yields the estimate 

Pg {x)<\x\- 2q '^ 
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with the exponent s := — 2q/ (2q — l) being such that 



s<-3/2 for l<g<3/2, s > -3/2 for q > 3/2. (9.1) 

We split the estimate for / by choosing g = l + l/A;>3/2 for \x\ < 1 and q e 
[1,3/2] for \x\>l so that 




\x\ 


< 1 A 


u| 


<V 1+lfk (\x\) 


\x\ 


< 1 A 


u| 


>V 1+1 f k (\x\) 


\x\ 


> 1 A 


u| 


<vM), 


\x\ 


> 1 A 


u| 


>V q (\x\). 



f(x,v)<g(x,v) : = C 
By (O, 

J Cr Sl with si > -8/5 for \x\ < 1, 
P f ~{x) < p g (x) < | CrS2 w . th S2 < _8/5 for \x\ > 1. 

By the Hardy-Littlewood-Sobolev inequality this implies that g has finite 
potential energy. The crucial step in the existence proof for a minimizer is 
the convergence of the potential energy along a minimizing sequence (fj). 
Since 0<fj< g, the finiteness of the potential energy for g allows us to pass 
to the limit using the dominated convergence theorem. 
Saturation of the constraints. 

Minimizers of the decoupled problems always saturate the constraints, i.e., 
||/o D ||i = ^; I |/o D | \i+i/k = N, and similarly for the flat case, because if for a 
minimizer one (or both) equalities were replaced by strict inequalities, then 
this minimizer can be rescaled in such a way that the constraints become 
saturated but the energy strictly decreases, which is a contradiction. A 
similar argument was used in the proof of Proposition 17.11 
The Euler- Lagrange relation and symmetry. 

For minimizers of the flat or non-flat problem the phase space distributions 
are functions of the local energy, more precisely, they satisfy relations as 
stated in Theorem I7.2[ the only differences being that 



E(x,v):=^\v\ 2 + U™(x), E{i,v):= l -\v\ 2 + U^ h {x) 



and the interaction term J Uq D Pq L in the definitions of Eq and Eq is dropped. 

The asserted symmetry of the minimizers follows from the fact that sym- 
metric decreasing rearrangements in x or x strictly decrease the energy except 
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when /9q D and Pq L and hence also the induced potentials are symmetric with 
respect to some point, cf. [151 3.7, 3.9]. 

Virial identity and compact support. 

Both flat and non-flat minimizers satisfy the virial identity 

2-Ekin(/ f f L ) =E pot (fQ L ), 2E kin (/o D ) = E pot (/o D )- 

This follows from the fact that these minimizers are time-independent solu- 
tions of the corresponding Vlasov-Poisson systems. A direct proof based on 
their minimizing property and scaling is given in [21 3.2]. The virial identities 
together with the restrictions on k and k immediately imply that the cut-off 
energies E and E Q in the Euler-Lagrange relations are strictly negative. In 
order to show that the minimizers have compact support it therefore suffices 
to show that their potentials converge to zero at spatial infinity. We indicate 
the corresponding argument for the 3D case, the flat one being completely 
analogous. Applying the Holder inequality to the first term in the estimate 

-CW</ r>o, 

J\x~y\<R \ x ~y\ ft 

implies that the potential is bounded and vanishes at spatial infinity, provided 
Pq D gL p (IR 3 ) with p>3/2. While a-priori this need not be the case for 0< 
k<7/2 we can use the fact that similar to (17. 3p . 

pf = c(E -U™)l, 

start with the known integrability f/ 3D GL 6 (l 3 ) to conclude that Pg D G 
L 6 / n (M 3 ), and obtain a new integrability estimate for the potential through 
the weak Young inequality. After finitely many iterations the desired inte- 
grability of pl B follows, cf. [22, Prop. 2.7]. 

Uniqueness for the 3D problem. 

First we notice that by the virial relation the Lagrange multipliers are 
uniquely determined by the constraint parameters M and N. In the 3D 
case we can now continue as follows. A minimizer is completely determined 
by its potential. The latter satisfies the Emden-Fowler equation 



l(r 2 U'y = c(E -U)l, i.e., ^y')' = -cy n + 



where y = E — U, and the constant c depends only on k and M, N. The 
solutions of this equation which are regular at the origin are uniquely de- 
termined by their value at the origin. Moreover, the scaling y a (r) =ay(a x r) 
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with A = (/c + l/2)/2 turns solutions into solutions. But the mass constraint 
fixes this scaling, and uniqueness of the minimizer follows. For more de- 
tails we refer to [221 p. 464]. Unfortunately, there is no analogue to the 
Emden-Fowler equation in the flat case — the flat potential does not satisfy 
the Poisson equation — and uniqueness in the flat case is not known. 
The radius relation \4-l\j in the 3D case. 

Each minimizer is a spherically symmetric steady state (f,p,U) of the three 
dimensional Vlasov-Poisson system, and for each choice of M and N there 
is a unique such steady state. If (f,p,U) is a steady state and a,/5>0 are 
arbitrary, then 

f a p(x,v ) = a 2 Pf(ax,(3v ), p a p(x) = ct 2 f3~ 2 p(ax) 1 U a p(x) = (3~ 2 U(ax) 

defines another one. There is a unique steady state (f*,p*,U*) with ||/*||i = 
1 = 1 1/*| |i+i/fe, and the minimizer with general M and N is obtained by rescal- 
ing /* with the parameters 

a = M (1 ~ 2fc)/3 iV 2fc+2)/3 (3 = M (fc ~ 2)/3 j\H fc+1)/3 . 

Since R = R*/a this implies (14.11) . 

The radius relation ^ n the flat case. 

We do not know whether for each choice of M and N there exists a unique 
minimizer /^ L , and so we cannot use the argument above to prove (14.21) . 

To obtain a two-parameter family of minimizers which obeys the radius 
relation (14. 2ft we proceed as follows. Since minimizers a-posteriori saturate 
the constraints we redefine 

Fm,n := { / G ^ (M 4 ) 1 1 1 /| | ! = M, 1 1 /| 1 1+1/S = N, E kin (f) < oo } . 
For / G J-~f\ we define the rescaled function 

fi*,u(x,v) : = pf(px,isv). 

Then 

||/^||i = M-^- 2 , ll/^IIi+i / S = M a - £)/Cl+i) ^- 2S/C * +1) , W(/^)=A*-^- 
For M, N > arbitrary we choose p, v such that f^ u 6 J-'f^, i.e., 

p = M-~ k N k+ \ v = M^ -1)/2 iV~^ +1)/2 . 
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The mapping !Fj\ 3 f i— > G .Fjy - N is one-to-one and onto. Let / denote 
an arbitrary, spherically symmetric minimizer of Tt over the set Then 
/J^ is a minimizer of 7i over JFj^y which we denote by ff^ N . To see this we 

observe that any function g G ^j^at can be written as <? = f^ u where / G J-fi, 
and hence 

n{g) = n(f^) = ^v-'nCf) > ^- 4 n(f h ) = 

If R* denotes the radius of the spatial support of / FL , then the spatial support 
of ff^ N has radius 

R = li - 1 R* = R*M-~ k N k+ \ 
and this is the remaining assertion (I4.2p . 
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